Theoretical Aspects of Gravitational Lensing in TeVe S 
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Since Bekenstein's (2004) creation of his Tensor- Vector-Scalar theory (TeVeS), the Modified 
Newtonian dynamics (MOND) paradigm has been redeemed from the embarrassment of lack- 
ing a relativistic version. One primary success of TeVe S is that it provides an enhancement of 
gravitational lensing, which could not be achieved by other MONDian theories. Following Beken- 
stein's work, we investigate the phenomena of gravitational lensing including deflection angles, 
lens equations and time delay. We find that the deflection angle would maintain its value while 
the distance of closest approach vary in the MOND regime. We also use the deflection angle law 
to derive magniflcation and investigate microlensing light curves. We find that the difference in 
the magnification of the two images in the point mass model is not a constant such as in GR. Be- 
sides, microlensing light curves could deviate significantly from GR in the deep MOND regime. 
Furthermore, the scalar field, which is introduced to enhance the deflection angle in TeVeS, 
contributes a negative effect on the potential time delay. Unfortunately this phenomenon is un- 
measurable in lensing systems where we can only observe the time delay between two images 
for a given source. However, this measurable time delay offers another constraint on the mass 
ratio of the dark matter and MOND scenarios, which in general differs from that given by the 
deflection angle. In other words, for a lensing system, if two masses, rrigN and mgM^ are mutually 
alternatives for the deflection angles in their own paradigm, regarding the time delay they are in 
general in an exclusive relation. 

Subject headings: gravitational Lensing — MOND — dark matter — gravitation — relativity 
1. Introduction velocities of individual galaxies are so large that 



ABSTRACT 



While Newtonian gravitational theory was ap- 
plied to the extragalactic region, it no longer 
passed the trials as easily as it had done in so- 
lar system. Within Coma and Virgo cluster the 



its total mass must exceed the sum of masses of 
its shining galaxies if the clusters are gravitational 
bound systems (Zwicky 1933; Smith 1936). The 
21 cm emission line from HI of spiral galaxies also 
show the need of extra masses besides visible stars 
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(van Albada ct al. 1985; Begman 1989). More- 
over, X-ray observations and strong gravitational 
lensing in early type galaxies and clusters imply 
that there are 60 ^ 85% masses in some unde- 
tected form (Fabbiano et al. 1989; Bohringer 1995; 
Loewenstein & White 1999). This so-called "miss- 
ing mass" problem may be explained in two dif- 
ferent ways: the existence of Dark Matter, and 
a modification of Newton law , hence of General 
Relativity (GR). The most famous representative 
of the modified gravity theory party is Modified 
Newtonian Dynamics (MOND) proposed by Mil- 
grom in 1983 (Milgrom 1983a, b,c), which assumes 
that Newton's second law needs to be modified 
when the acceleration is very small: 

/i(|a|/ao)a= -V$N- (1) 

Here $n is the usual Newtonian potential of 
the visible matter and Oq ~ 1 x 10~^° m s~^ 
from the empirical data, such as the TuUy-Fishcr 
relation and rotation curves (Sanders & Verheijen 
1998); tl{x) « a; for a; -C 1, fi{x) ^ 1 for x > 1. 
In the solar system where accelerations arc strong 
compared to ao, the formula (1) is just Newton's 
second law a = — V$n- 

From this simple formulation the TuUy-Fisher 
relation and flat rotation curves follow triv- 
ially (Sanders & McGaugh 2002). The rotation 
curves of spiral galaxies, especially of LSBs, can 
be explained extremely well (Begeman et al. 1991; 
Sanders 1996; McGaugh & de Blok 1998; Sanders 
& Verheijen 1998). Moreover, the dearth of dark 
matter in some ordinate elliptical galaxies can be 
intrinsically explained in MOND paradigm (Mil- 
grom & Sanders 2003). Even though many such 
phenomena in galaxy systems can be explained 
by this simple modification with an astonish- 
ing precision, the MOND paradigm had lacked 
a flawless relativistic gravitational theory for two 
decades (Bekenstein & Milgrom 1984; Bekenstein 
1988; Sanders 1988, 1997), and thus could not 
address the lensing and other relativistic gravi- 
tational observations. However, since the end of 
March 2004, the embarrassment has been removed 
by Bekenstein (2004), when he finally contrived 
a Lorentz covariant relativistic gravitational the- 
ory called TeVeS, which is based on three dy- 
namical gravitational fields: an Einstein metric 
g^v, a vector field il^ and a scalar field (p (the 
acronym TeVeS is for this Tensor- Vector-Scalar 



content). In this scheme, TeVeS contains MOND 
for nonrelativistic dynamics, agrees with the solar 
system measurements for the f3 and 7 PPN coeSi- 
cients and avoids superluminal propagation of the 
metric, scalar and vector waves. Moreover, while 
concerning cosmogony, it not only has cosmolog- 
ical models which arc similar to that of GR but 
also overcomes the Silk damping problem of the 
CDM-dearth universe with the help of the scalar 
field (Skordis et al. 2005). Now, we can investigate 
the relativistic MOND paradigm on many obser- 
vations and experiments (Giannios 2005; Skordis 
et al. 2005; Zhao 2005a,b; Zhao et al. 2005). 

On the other hand, the GR with Dark Matter 
(DM) paradigm, despite its success for the large 
scale structure and the cosmic microwave back- 
ground (CMB)(see e.g. Spergel et al. 2003), has 
some problems on the galactic scales. For ex- 
ample, the dark matter halos, if they really ex- 
ist, must have a shallower core than the univer- 
sal NFW density profile had predicted (cuspy core 
problem, see e.g. Navarro & Steinmetz 2000; Ro- 
manowsky et al. 2003); the disc size will be too 
small while considering the angular transforma- 
tion from baryon to dark matter halo (angular 
momentum catastrophe, see e.g. Navarro et al. 
1995; Moore et al. 1999; Mailer fc Dekel 2002; 
Kazantzidis et al. 2004); the number of sub-halos 
predicted by non-baryon simulation is much more 
than the observed number of dwarf galaxies (sub- 
structure problems, see e.g. Kauffmann et al. 1993; 
Navarro et al. 1995). Despite these problems, the 
GR with DM paradigm is less falsifiable and sup- 
ported by a larger group because of its unlimited 
choice of models and parameters. 

In order to judge which side better accounts for 
the phenomena of the natural world, the dispar- 
ity between TeVe S and GR with DM paradigms 
has to be closely investigated. This distinction 
will be important for the physics of MOND and 
the "missing mass" problem. Gravitational lens- 
ing (GL) phenomena provide us an opportunity 
because they dominate on scales that even surpass 
the size of dark matter halo (Kochanck 2004). 

For a long time, light bending has been a great 
challenge for Modified Gravity. Any scalar field 
jointly coupling with the Einstein metric via a con- 
formal transformation would yield the same opti- 
cal phenomena as in the standard Einstein metric. 
Thus in the context of MOND, where little dark 
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matter exists, it is hard to produce the anoma- 
lously strong deflection of photons (Bekenstein & 
Sanders 1994; Soussa fc Woodard 2003, 2004). 
However, this problem can be solved by postu- 
lating the disformal relation between the physi- 
cal metric and the Einstein metric (Sanders 1997). 
This is what Sanders' stratified theory, and then 
TeVe S take over. 

In this paper, we will discuss the theoretical as- 
pects of Gravitational lensing systems within the 
framework of TeVeS. In § 2, we summarize the 
characteristics of TeVe S and introduce a cosmol- 
ogy model within which gravitational lensing phe- 
nomena will be discussed. In § 3 we derive the 
deflection angle in the general form and its two 
particular limits: the Newtonian and the MOND 
regimes. From this we obtain lens equations, 
which are the foundation of the other GL phe- 
nomena. In § 4, we apply the standard time delay 
formalism in TeVe S. Finally, § 5 summarizes the 
presuppositions of this paper and discusses some 
possible conclusions of our results in § 3 and § 4. 

2. Fundamentals of TeVe S and The Cos- 
mological Model 

Before introducing the action of TeVeS, we 

should explain some important features of the dis- 
formal metric relation. In Bekenstein's Tensor- 
Vector-Scalar theory (TeVeS), the role of Ein- 
stein's metric gap is replaced by the physical met- 
ric (Bekenstein 2004) 

~9a0 = e-^*{ga0 + ila%) - e'^'^Haiiff, (2) 
with a normalization condition on the vector field: 

5"^ilail/3 = -1. (3) 

The physical metric, which governs the dynamical 
behavior of matter (the equivalence principle i.e. 
minimal coupling ) is simultaneously influenced by 
a 4- vector field ila, a scalar field cp and the Ein- 
stein's gravity gajs- 

From the theoretical viewpoint, if DM does not 
exist in galaxies or cluster scale, the concept of 
geodesies must be different from that of GR. Al- 
though there is no deeper theoretical base for the 
disformal relation (2), empirical testament like 
lensing effects will be important to justify this 
choice (Sanders 1997). 



The cosmological model, which determines the 
background spacetime structure, is an essential in- 
gredient in gravitational lensing (Schneider et al. 
1992; Fetters, Levine & Wambsganss 2001). On 
the one hand, the distance in the gravitational 
lensing scenario would vary under different cosmo- 
logical models. On the other hand, gravitational 
lensing provides a tool to determinate Hq and to 
constrain the cosmological parameters (Blandford 
& Narayan 1986; Schneider 1996). Since we have 
used an alternative gravitation theory, it is impor- 
tant to discuss its effects on cosmology, which is 
solely influenced by gravity. In this section, we in- 
troduce Bekenstein's Relativistic MOND theory, 
and then describe the revised Friedman model. 

2.1. Actions of TeVe S 

The convention in this paper are the metric sig- 
nature +2 and units with c = 1. 

There are three dynamical fields gaff, and (p 
built in TeVe S. In addition to these fields, there 
is another non-dynamical scalar field a, which is 
introduced into the scalar fleld action to eliminate 
the superluminal problem of (Bekenstein 2004). 

By varying the actions for these fields, the re- 
spective equations of motion for the scalar fields, 
vector fields and the Einstein metric can be arrived 
at, therefore the physical metric is determined. 

The action in TeVe S can be divided into four 
parts, S = Sg + Ss + Sy + Sm- The geometrical 
action, 

Sg = {16nG)-' J g'^^R^0{-g)'/^d^x, (4) 

is identical to the Einstein-Hilbert action in GR 
and is used to produce the Einstein metric, gaf}- 
The matter action coupling to the scalar field, (j), 
and vector field, iia, through (jap is taken to be 

Sm = J C{-~gY'^(fx, (5) 

where £ is a lagrangian density for the fields under 
consideration and should be considered as func- 
tional of the physical metric and its derivatives. 

The vector action with the form 

= "3^ / 

-2{X/K){g'^^^^'d, + 1)] (-ff)V2rf4^6) 
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yields the vector field equation. Here, X is a di- 
mensionless parameter, and A is a Lagrange mul- 
tiplier. Moreover, il[a,/i] means iia,^i — ^^i,a- 

The action of the two scalar fields is taken to 
have the form 



\a<P,l3 



+ ^Gr^a*FikGa^)]i-g)'/^d*x, (7) 

where h"'^ = g"'^ — if"!!'': F is a free function 
in order to produce the dynamical behaviors of 
MOND, which needs to be constrained empirically 
from cosmological model; I is a constant length 
and k is another dimensionless parameter in this 
theory. In the FRW-like cosmological models, k is 
smaller than 10~^ (Bckenstein 2004; Skordis et al. 
2005). Variation of a of the action gives 



-,,F{^,)-^|,'F'{|,)=y 



(8) 



where y = kf.'^h'^'^(f>,f^(t)^,y is a free fimction of fi 
equivalent to F. However, in later discussions, it 
is more convenient to consider fi{y) as a function 
of y. Variation of ^ of the action then gives 

+ (1 + e-^^)il«il'5]T„;3. (9) 

As a through n = kGa'^ can be expressed in terms 
of (j) explicitly, Eq. (9) involves ^ only. 

To illustrate ideas, in this work we adopt the 
free function suggested by Bckenstein (2004) 



y=4 



3/i2(/i-2)2 



1 



(10) 



whence y approaches zero asymptotically as y ~ 
when — > 0, and y — > oo when — > 1. 
We are going to discuss gravitational Icnsing of 
a point mass in § 3 and § 4. Thus we are interested 
in static and spherically symmetric geometry. In 
this case ^ runs from 1 (the Newtonian regime) 
to (the deep MOND regime). In order to sat- 
isfy the behavior of the potential at these limits 
we require y gc asymptotically as ^ 0, and 
y — > (X) as ^ 1. In the intermediate regime j/(/ti) 
(or /i(y)) is arbitrary. Nonetheless, the exact func- 
tional form does not affect our results because we 
do not consider the intermediate regime. 



2.2. The Cosmological Model 

As pointed out by Bckenstein (2004), the phys- 
ical metric of FRW cosmology in TeVeS is ob- 
tained by replacing dt — » e'^dt and a — > e~'^a in 
the primitive Roberson- Walker metric. Here, the 
Friedmann equation becomes: 



a \a 



(11) 



where a/a and a/ a are essentially equal in the 
whole history because e^*^ ~ 1 and ~ (Bcken- 
stein 2004; Hao & Akhoury 2005). 

Since we have the additional scalar fields and 
vector fields in Bekenstein's picture, the Einstein 
equation, which tells us the geometrical prop- 
erties, is not only determined by the energy- 
momentum tensor as before, but also by the scalar 
and vector fields. Here, <p, a, and il^ are assumed 
to partake of the symmetry of the spacetime. The 
vector field is considered to be parallel with the 
cosmological time, i.e. il" = 6", because there 
should be no preferred spatial direction (Bcken- 
stein 2004). Moreover, if we apply Weyl's postu- 
late, which states that galaxies would move along 
nonintersecting world lines and therefore the spa- 
tial components in the energy-momentum tensor 
can be neglected (Weyl 1923), then only the tt 
component of Einstein's equation remains: 



27r 



m2F(m)+A 
47r 



lyF{n)+'-^^'F'{n) 



+ A. (12) 



The appearance of A in Eq. (12) results from the 
choice of the free function F(p), which is decided 
by Eq. (8) and Eq. (10). It is also possible to con- 
struct another form of F{n) that makes the dark 
energy as an effective phenomenon of the scalar 
field (j) (Hao & Akhoury 2005). However, before 
having any theoretical reason for the choice of the 
free function F{fi), we would like to leave the in- 
terpretation alone and only treat it as a parameter 
in the cosmological model. 

By the Roberson- Walker metric, which deter- 
mines the tt component of the Einstein tensor 
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Eq. (12), we can get the modified Friedmann equa- 
tion in a flat universe in the form 



a? SttG , SttG 



2 — Pm -\ —P<l> + 



(13) 



where pM = pe corresponds to the matter den- 
sity in Einstein's picture and 



P4> 



[^f^'F{„) + 'yF'{„)] (14) 



is the density due to the scalar field. For con- 
venience, we can rewrite Eq. (13) in terms of Q 
density: 



(15) 



Since gravitational lensing must be observed 
within the framework of galaxies or clusters of 
galaxies, we only need to consider the cosmologi- 
cal model in the matter and A-dominated era, so 
Slm{z) — f^m(l + z)^- Furthermore, in order to 
avoid a significant effect of the integrated Sachs 
Wolfe term on CMB power spectrum, n^{z) must 
be extremely small compared with the matter 
density (Skordis et al. 2005) . So even if we do not 
know exactly how n^{z) evolves with z, we can 
neglect it, because its value is overwhelmed by the 
uncertainty of Q.^ and Q\ for a non-CDM uni- 
verse (McGaugh 1999, 2004; Skordis et al. 2005). 

To summarize, in TeVe S the angular distance 
at any redshift can be calculated by 

dA ^ ^-^-^TT^ r dz p^l + zf + n^] 
-no Jo 

(16) 

We conclude that if there is any difference in the 
angular distance between GR and TeVe S, it arises 
from the low matter density universe in TeVe S 
rather than the influence of the scalar field ^. 

Finally, it is worth noting that although it is not 
trivial to consider an inhomogeneous universe on 
the gravitational lens scale (Dyer & Roder 1973), 
we do not consider that case because wc only deal 
with a single point mass lens embedded in a back- 
ground Friedmann universe in this paper. 

3. Gravitational Lensing in TeVe S 

When building gravitational lensing models, 
people usually use some approximations, such as 
the static and thin lens approximations. We 



also adopt the weak field assumption since all 
known gravitational lens phenomena are influ- 
enced by weak gravitational systems. Although 
the Schwarzschild lens (point mass model) is not 
sufflcient for gravitational lensing effects, we con- 
sider this for its simplicity. Moreover it gives re- 
sults of the same ordcr-of-magnitudc as those for 
more realistic lenses (Schneider et al. 1992). 

3.1. Deflection Angle in Terms of Mass 

In the weak field limit, the physical metric of a 
static spherically symmetric system in TeVeS can 
be expressed in the isotropic form: 

5c/3 dx°'dx^ = -(1 2^)dt^ 

where $ = S$Ar + </) and S = e-2'^<=(l KI2)-^ 
with (\)(. as the asymptotic boundary value of (/> and 
K < 10~^ from the constrain of PPN parameters 
(Bekenstein 2004). 

Considering a quasistatic system with a perfect 
fluid, the scalar equation (9) can be reduced into 



V- M(y)V?!> =kGp 



(18) 



because iia has only a time-independent tempo- 
ral component and ^ ^ 1 in this case. Here 
y = fc£^(V0)^, which is obtained from the general 
form, y = k£'^h^'^(j)^^(j)^^, in spherical symmetric 
geometry. Prom Eq. (18), we can find a relation 
between $jv and ^, 



V0 = (k/Ann) V$jv 



(19) 



Therefore, $ can be determined whenever /x is 
known. 

It is also worth to stress that gu = — (1 + 2$) 
follows the empirical definition of the gravitational 
mass nig in the solar system when fi ^ 1: 



9 = 



1 dgtt 

2 dg 



(S + fc/4.)^ = ^ 



(20) 



and that the spatial components, gu = (1 — 2$), 
guarantee the theory passes the classical post- 
Newtonian (PPN) tests of gravity (Bekenstein 
2004). 

We now apply the static spherically symmetric 
metric, c.f. Eq. (17), to the GL systems. Con- 
sidering a light ray, which propagates on a null 
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geodesic, and moves in the equatorial plane, i.e. 
= 7r/2, we can arrive at three constants of mo- 
tion from the physical metric in equation (17) with 
some proper affine parameter, r: 

(l + 2$)i = £;, (21) 

(l-2$)eV = i, (22) 
(1 - 2$)^2 + (1 + 2^)g-'^L^ - (1 - 2$)i;2 = 0, (23) 

where over dot denotes the derivative with re- 
spect to T. Eliminating E and L in Eq. (23) from 
Eq. (21) and (22), and recalling the fact that g 
vanishes at the closest approach, g = go, yields 



L2 
£2 



go(l-2$o) 
(l + 2$o) ' 



(24) 



where, $0 = ^(qo)- Combining Eq. (24) and 
the three constants of motion, gives the equation 
which describes the shape of the photon orbit: 

-il~A^) + il-A^o)igo/gf[g-^{dg/dipf + l] = 0. 

(25) 

The solution of this equation in quadrature is 



r|(^)Mi-4(<i>-$o)]-ij"''^ 



(26) 

If we take the Taylor expansion of Eq. (26) to 
the first order of $, then it is easy to see that 
the quadrature is a combination of the angle of 
a straight line feeling no gravity and the angle of 
the deviation due to the gravitation field. Hence 
the deflection angle can be approximated from the 
first order terms of this Taylor expansion: 



A(/3 = go 



2£)($ - $0) 



£>o {g^ - £»5) 



dg. 



(27) 



Furthermore, after some manipulations (see e.g. 
Bekenstein 2004), the deflection angle is 



= go 



4$' 



dg - go- 



4$ 



[g2-gl)i/2 - ^"{g^~glY/^- 

(28) 

It is obvious that the deflection angle between 
the closest approach and any position is always 
positive. First, gravity is always attractive (there- 
fore > 0), so the value of the quadrature is pos- 
itive because every term in it is positive. Second, 
the second term of the right hand side is negligible 
even in the MOND regime, in which $ behaves as 
In g, if g is large. 



MOND Regime 




intermediate 

MOND 
1^ Regime 



-3 



-2 



bq p /r 
^ in ' n' n 



Fig. 1. — Deflection angle of light derived under 
the point mass model in TeVeS with k = 0.01 
and K = 0.01(solid line) and GR (dashed line). 
The closest approach is normalized to the length 
scale , which depends on the mass of the lens. 
The dotted line shows the work of Mortlock and 
Turner (2001), who started from a reliable intu- 
ition without a relativistic MOND. In the two ex- 
treme limits, their result is exactly the same as the 
deflection law derived in TeVeS. However, in the 
intermediate region there is a difference. 

3.2. Two Limits of the Deflection Angle 

Even though we can make sure that the de- 
flection angle can be enhanced in TeVe S by the 
simple argument above, we need the exact form of 
the the potential to reach the real value of a de- 
flection angle. Unfortunately, the form of poten- 
tial depends on the free function F(fi), which can 
only be constrained empirically. However, there 
are two limits, the Newtonian and the MOND 
regime, in which the dynamical behavior has been 
clearly known, and so to the potential form. Ac- 
tually, any choice of the F{fi) cannot contradict 
the empirical potential form in these two limits. 
Here, we discuss the deflection angles in the two 
regimes (see Fig. 1). 

The MOND regime means V0 is significant 
compared with V$, whence by Eq. (19), /it <C 1 
has to be a sufficient condition for this regime. 
From Eq. (10) and y — fc£2(V0)2, we can obtain 
^l « (fc/3)i/2Z|V0| when n < 1 (i.e., in MOND 
regime). It has to be mentioned that although 
this relation results from the freedom of an arbi- 
trary function, /i ~ |V0| should always persist in 
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order to guarantee Vi/) varying as 1/ g m the deep 
MOND regime, a condition which follows from the 
MOND paradigm. Using this relation in Eq. (19) 
to eliminate V(/> and replacing V$jv with V$ gives 



M = (A:/87r5)(-l+v'l + 4|$'|/aoj = — I — 
Here, the constant 



(3A:) 



1/2 



A-Km 



(29) 



(30) 



can be identified as Milgrom's constant. We 

should keep in mind that this form of /i is valid 
under the condition <g; (47r/fc)^ao. It also offers 
a criterion for distinguishing the Newtonian and 
the MOND regime. 
We define a distance, 



G 



Nnig 



1/2 



(31) 



at which acceleration, a, equals to Milgrom's con- 
stant, Oq. Then a particle is in the deep MOND 
regime if :2> kro/An, and in the Newtonian 
regime if ^ -C kro/Aw. 

3.2.1. The Newtonian regime 

Although a photon moving in a GL system may 
come from a distance q ^ kro/An to the closest 
approach Qq kro/Air and then fly away again, 
the influence of gravity only dominates in the re- 
gion which is very near the gravitation center. 
Therefore we call a GL system as being in the 
Newtonian regime if ^ kro/Air. 

In the Newtonian regime, in Eq (28) can 
be calculated from the relation i> = S$ n + 4' ^iid 
Eq. (19) for — > 1. This would yield the deflection 
angle of a Schwarzschild lens. 



k ^ AGNmg 
Att go 



(32) 



Here again, Gjv is identical to the gravitational 
constant measured in a local experiment. This is 
the same as the result from GR. 

3.2.2. The Deep MOND Regime 

In the deep MOND regime, go ^ kro/Air, n has 
the form of Eq. (29). Then along with Eq. (19) 



Fig. 2. — Embedding diagram of light bending. 

For a thin lens, the lensing equation and geomet- 
rical time delay can be immediately read off. 

and Eq. (29), the deflection angle for a point mass 
model in the deep MOND regime can be arrived 
at from Eq. (28): 



go{^ + k/Air) 

"ao(S2 + fcS/47r) 




GNrUg 



1/2' 



,(33) 



with S = 1 — K/2 — 20c- Since all of K, (f)^ and 
k are much less than 1, for a given mass, the de- 
flection angle in the deep MOND regime differs 
from that in GR (or equivalence in the Newtonian 
regime; Eq. [32]) by an amount almost indepen- 
dent of the distance of the closest approach. 

We have to stress that even though we do not 
know exactly how deflection angle varies with re- 
spect to the closest approach in the intermedi- 
ate MOND regime, our result based on the deep 
MOND assumption (i.e., fi <C 1), which is incor- 
rect in the regime go < kro/Air, approaches to 
the Newtonian prediction while decreasing (see 
Fig. 1). Therefore we believe even in the interme- 
diate regime, the deflection law based on the deep 
MOND assumption is approximately correct. 

3.3. Magnification and Microlensing 

From the observational point of view, the mea- 
surable data is not deflection angles but posi- 
tions of the projected sources, i.e. 6 = go/Dj, . 
Here, D/, is the angular distance of lens (Eq. [16]). 
Therefore we need a relation to connect Q and Ai^, 
the deflection angle derived in the last section. We 
can obtain their relation from an embedding dia- 



7 



gram (Fig. 2) and get the lens equations: 

P = e^-a{e^), (34) 
l3 = a{e_)-0_. (35) 

Here the reduced deflection angle is defined as 

a{e) = {Dls/Ds)Aip{9). In GR it is not impor- 
tant to distinguish the difference between Eq. (34) 
and Eq. (35) in finding 9 Paczyhski (1986). It is, 
however, not the case in TeVe S. In order to keep 
0+ = d_ when /? ^ 0, we must use Eq. (34) for 
and Eq. (35) for 9_ . 

After substituting A(p and 9 into the lens equa- 
tion Eq. (34), wo get 



P = 



Dls 4:GNmg 
DlDs 9 

: ttODl 



E + k/An 2 ro ^ An' ^ 



(36) 



2 



Fig. 3. — The difference in the magnification of 
the two projected images. In TeVe S, it deviates 
from the constant vahic in GR (dashed line) . From 
the solid to dottcd-dashed line, they represent the 
cases of TeVeS with 9e/9o = 1.0, 0.5 and 0.1. 



For clarity, we can define 6q = ro/DL and 6e with 
the form 



E = AGNmg 



Dls \ 
DlDsJ 



1/2 



(37) 



9e is called the Einstein radius; it is used as a 
length scale in the standard GL modeling. As 
fc <C 1 and <C 1, the lens equation (36) can 
be reduced to 



6 



■[1 + 



TT 9 . 

2^^ 



with the reduced deflection angle, a{6), as 



2 9a 



(38) 



(39) 



If Eq. (35) is used instead of Eq. (35), then the 
only difference in Eq. (38) is from /3 to — /3. It is 
easily seen that the effective contribution of the 
scalar field can be expressed as 1x9 129^ and is 
independent of 9q (or equivalently ^o)- 

For a point source, there are always two 
projected images, which can be obtained from 
Eq. (34) and Eq. (35). The solutions are 



^± = ±7 



1 



^ 29a 



± 



+ 46»' 



1/2^ 



(40) 



Moreover, since gravitational lensing conserves 
surface brightness, the magniflcations of the image 
to source intensity are determined by their area 
ratio, 

9± d9± 



d0 



(41) 



Interesting enough, the difference in the magnifi- 
cation of the two images is no longer constant in 
TeVe S as in GR (see Fig. 3): 



> 1 



(42) 



And in the deep MOND regime, the total magni- 
fication is given by 



Ui (ui2 +4) 



1/2 



+ 



\2uJ 



U2 



U2 (^2^ + 4) 



1/2 



where 



/? n 9e n 9e 

U = — ; Ui =«+-— ; U2 = U-- — . 



(43) 
,(44) 

(45) 



The magnification reduces to the standard GR 

form as 9o (or tq) — > oo. It is then straight- 
> Oforward to study the light curves of microlens- 
ing (Paczyhski 1986). 
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Fig. 4. — upper panel: Microlensing light curves 
in TeVeS for 6'is/0o = 1-0. Prom solid line to 
dashed line, fimin/f^E = 0.1, 0.3, 0.4, 0.5. lower 
panel: Light curves of microlensing in TeVeS 
with Oe/Oq = 0.1 (dashed-dotted), 0.5 (dotted), 
1.0 (solid) and in GR (dashed line). The differ- 
ence of the magnitude is given by Am = 2.5 log A 
and to is the time that sources pass through Ein- 
stein radius. 



Microlensing events results from the relative 
proper motions of lenses to the sources. Be- 
cause of proper motions, the alignment between 
lens-observer and source-observer would vary with 
time, 



u{t) = ^/iP~j9ErTWtoy 



(46) 



and so does the magnification. Here (3min oc- 
curs at the greatest of the alignment and to is 
the time it takes the source to move with re- 
spect to the lens by one 9e- Putting Eq. (46) 
into Eq. (44) and Eq. (45), we can obtain the mi- 
crolensing light curves expressed in stellar magni- 
tudes (Am = 2.5 log A). Like the case in GR, the 
variation reaches the maximum as f5{t) = fimin, 
and for any given 7r0g^/20„, the smaller the (3minj 
the greater amplitude of the light curve. However, 
the shape of light curves in TeVe S would deviate 
from that in GR as Tr6^^/26g is significant (i.e., in 



deep MOND regime). In general the contribution 
of 7r6»gV26'o would raise the peak (see Fig. 4). 

Even though the differences in microlensing 
events between TeVeS and GR only occurs at 
source redshift Zg > 1 (Mortlock & Turner 2001), 
which is beyond the current micolensing projects 
within the local group (see e.g., Alcock et al. 2000; 
Paczyhski 1996), observations such as several 
quasar microlensing events (see e.g., Wambsganss 
2001) may help us to decide whether MOND is a 
viable alternative to the dark matter paradigm. 

4. Time Delay 

Time-delay always plays an important role in 

gravitational lensing. It not only provides a more 
convenient way to get the lens equation (Schnei- 
der 1985; Blandford & Narayan 1986) but also 
has some important values in its own right (for 
more complete reviews, see Courbin et al. 2002 or 
Kochanek & Schechter 2003). As it is an alter- 
ative to the cold dark matter, it should be very 
interesting to study the behavior of time delay in 
TeVeS. 

4.1. Arrival Time in TeVeS 

We apply the static and isotropic metric de- 
scribed by Eq. (17) and assume that the photons 
are confined to the (equatorial planc^ without loss 
of generality . Then, following the standard pro- 
cedure of light bending (see e.g. Weinberg 1972), 
not only the deflection angle can be derived, but 
also the equation which describes the arrival time 
of a light ray 

This equation has a quadrature as a solution, 
which states that a photon from qq to any distance 
g (or from g to the closest approach distance go) 
would require a time, 



J 00 



(1 -4$)^/^de 



[l-(l+4$)/(l + 4$o)(f )2]i/2 



• (48) 



Taylor expanding Eq. (48) to the first order of 
0($), the original quadrature can be reduced to 



t{g 



,go)= / dg[l-{—) ] 
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+ 



/ 



s , 2($ 
dg 

eo 



*o)(eo/£')^ 



[1 



-/■ 



dg- 



(^)2]3/2 

2$ 



(M)2]l/2 



(49) 



The second quadrature of Eq. (49) differs from 
Eq. (27) only by a factor go, whence, from this 
part, the arrival time contributions due to both 
gravitational mass and the scalar field are posi- 
tive. However, this is not the case in the third 
quadrature of Eq. (49). Since $7v and share an 
opposite sign in their contributions to the resul- 
tant arrival time are also opposite. Moreover since 
the third quadrature is always larger than the sec- 
ond one (it is easily to observe this by multiplying 
by a factor "1 — {go/ q)^" in both numerator and 
denominator of the second quadrature), the time 
delay due to the scalar field will always be opposite 
to that due to the gravitational mass, i.e. with an 
opposite sign. In other words, the arrival time will 
be shorten rather than enhanced by the positive 
scalar field; this is a necessary condition in TeVe S 
for not violating causality (Bekenstein 2004). 

4.2. Arrival Time and the MOND Regime 

As in the case of the deflection angle, in or- 
der to calculate time delay from Eq. (49), we need 
some additional information on the free function, 
F{^i) to determine the potential i> in TeVe S. How- 
ever, there are some differences in these two con- 
siderations. The first one arises from the fact that 
time delay is influenced by gravity even far away 
from the center. This makes a purely Newtonian 
regime unavailable. Therefore, we only consider 
the MOND regime here. Moreover, unlike the de- 
flection angle case, we need information about 
in addition to for the time delay. Eliminating 
II of Eq. (19) from Eq. (29) yields 

<^ = (SaoGm3)i/2ln^ + ,^e, (50) 

where (/>c can be absorbed into the symmetric 
physical metric g^^ by rescaling the t and g coor- 
dinates appropriately (Bekenstein 2004). Recall- 
ing that <I> = S$jv + 0, and putting Eq. (49) into 
Eq. (50), we get the arrival time in the MOND 
regime: 



t{g,go) = {g'-gir' + 2^GmJ'—f) 

\Q+ QoJ 



1/2 



5x10'° 



-5x1010 



-1x10" 



0.0 



0.5 



1.0 
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Fig. 5. — The time delay (the real light path mi- 
nus the straight between sources and observers) 
of a Icnsing system with a Af = 5 x IQ'^'^Mq lens 
in the MOND regime. The closest approach is 
normalized to the length scale r^, and both the ob- 
server and the source are located at 10 r^. The 
dotted line and the dashed line show the potential 
time delay induced by the ordinary matter and 
scalar field respectively. Moreover the geometrical 
parts are represented by the dotted-dashed (scalar 
field) and dotted-dotted-dashed line( visible mat- 
ter). Adding all components yields the total time 
delay (solid line) and shows a negative effect on 
the arrival time in TeVe S. 



+2SGm,ln(VZZl±^) 



£•0 



-2 {EaoGmgf^ 



(ln(|^)-l)(,^-,^)V2 + ,osec-i(^) 
kro go 

+2 {EaoGmg 



a/2 



gosec-\^)-gllni^){g^-gl)-'/^ 



.(51) 



Here again, G differs by a factor of S -|- k/Aw from 
the usual Newtonian constant, which is measured 
in the solar system experiments. We would like to 
stress that as a particular case of our cinsideration 
in § 4.1, the sum of last two terms in Eq. (51) is 
always negative. 

Intrinsically, we may think that time delay 
(photon deviation of the path length) should in- 
crease with deflection angle as read off from the 
embedding diagram (Fig. 2). However, this is only 
true if we do not consider the deflection poten- 
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tial. Actually, a positive scalar field will offer a 
hyperbolic-like spacetime, in which the distance is 
shorter than that of a flat universe (see Fig. 5). 

4.3. Measurable Time Delay 

Since it is the positions of sources and lens, not 
the distance from sources (or observers) to lens 
that arc observed, time delay is conventionally ex- 
pressed in terms of these dimensionless angles. It 
is appropriate to rewrite Eq. (49) as: 



t 



I'-'-f 

•'Bo •'00 



2M, 



(52) 



here, dl = [1 - {1 + 4$/l + 4^o){f?]~^^'^dQ is a 
line element like that of Euclidean geometry, i.e. it 
obeys dl^ = dg^ + g^dif^ . Traditionally, we call the 
first quadrature in Eq. (52) the geometrical arrival 
time because it can be read off in a geometrical 
way (see Fig. 2): 



tgeom{Qi,Q2) = ^ [Qo - v)^ + Dig + ^ qI + (53) 

Here 77 = DsP. Indeed if wc let a photon move 
from §2 to go and from go to gi, it can be proved 
that the contributions to the arrival time cal- 
culated from the first two integrals of Eq. (49) 
are identical to the geometrical arrival time, c.f. 
Eq. (53). On the other hand, the third quadrature 
in Eq. (49), which equals the second integral of 
Eq. (52), contributes to what is called the potential 
time delay. Adding the geometrical and potential 
contributions to the arrival time and subtracting 
the arrival time for an unlensed ray from S to O, 
and considering the expansion of the cosmological 
background, wc obtain the time delay of a possible 
ray compares to an un-defiected ray. 



At={l + z) 



go 
Dl 



JL 

Ds 



)'-V'(go) 



-hi. 

(54) 

where, ip{go) called the deflection potential can be 
calculated from the potential time delay. Once 
again we only consider what happens in the 
MOND regime. For any given sources and lenses, 
the distances from the lenses to the observer (or 
sources) are decided, whence the potential time 
delay, which is given from the second and the fifth 
terms in Eq. (51), can be reduced to some constant 
plus the deflection potential contributions 
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Fig. 6. — The measurable time delay (time delay 
between two images for a given source) of a grav- 
itational lensing system in TeVe S. The system is 
assumed to have its largest closest approach at 
g^^ = 2rj,and a lens with M = 5 x 10^^ Af© to 
resemble a cluster of galaxies alike Virgo cluster. 
The dotted line is the contribution from the visible 
mass, and the dashed line expresses the contribu- 
tion from the scalar field. Their combination is 
showed by the solid fine. When bgj g^^ 0, the 
discrepancy between the time delay of the scalar 
field and ordinary matter decreases. On the other 
hand, when Sgj g^ — > 1, the infiuence of the or- 
dinary matter increases and ultimately exceeds 
that of the scalar field. For the same lens, two 
systems with g^ = 5 ''0 ^nd g^ = 10 respec- 
tively are presented by dotted-dashed and dotted- 
dotted-dashed line. 



aoGNrUg 



1/2 



4GNmg 



1 -I- 47rS/fc 



(56) 
(57) 



We have shown in § 4.2 that the time delay (rel- 
ative to the undeflected path) is reduced rather en- 
hanced in the TeVe S framework. Unfortunately, 
this phenomena cannot be measured in any GL 
system. For any given source at position /3, we are 
only able to measure the time difference for two 



images g'^^^ 



g\ ', which share the same lens po- 



V'GiJ = ^GNmAng^, 



(55) 



tential, lens-observer, and lens-sources distances. 
In other words, the constant B in Eq. (54) is the 
same for all light rays from the source plane to 
the observer. Therefore in regard to the potential 
time delay, we can only measure the contributions 
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Fig. 7. — Mass ratios of the MOND and dark 
matter paradigms. The soHd Hne shows the con- 
straint on the mass ratio while requiring an iden- 
tical amount of deflection angle in the two differ- 
ent paradigm. This constrain is found again from 
the requirement of the measurable time delay as 
5g ^ Qg. However, the corresponding mass ratios 
of the time delay will be increasingly less than that 
of the deflection angle as the difference between 
the two closest approach increases. From the dot- 
ted line to the dotted-dotted-dashed line, bI/bI = 
2 ,3, 4, 5. 



from the deflection potential, which possesses the 
same sign in GR and TeVeS (Fig. 6). Moreover, 
since the deflection potential does not depend on 
the length scale kr^/^n, the free choice of has 
no influence on the measurable time delay in GL 
systems. 

Even though we cannot measure the opposite 
contribution to the arrival time in TeVe S, time 
delay between two images due to the deflection 
potential offers another constraint on the mass ra- 
tio, which is not a priori identical to that due to 
deflection angle. 

In order to illustrate this, we consider a general 
case: the value of is "n" times g'^, whence the 
time delay between the two images is 



Aij; = AGNmg 



Inn -|- 



TT (n-1) 



(58) 



For the extreme case, the difference between g^^ 
and is negligible compared to either distance. 
Hence, ln(^i/£i^) ~ Sg^/ g^, and Eq. (58) can be 



5Bo 

' Bo 



1 



2 r„ 



(59) 



If we compare Eq. (59) with Eq. (39), we can find 
that the similarity of the two equations makes the 
mass ratio nigN / nigM exactly the same under the 
requirement of cither identical deflection angle or 
time delay. However, if the difference between g^ 
and varrho'^^ is not small, the mass ratio given by 
the time delay will be less than that obtained from 
the deflection angle (Fig. 7). 

4.4. Time Delay and Lensing Equation 

As showed by Schnieder (1985) and Blandford 
and Narayan (1986), time delay with Fermat's 
principle offers an alternative for deriving the lens 
mapping. Moreover, since Fermat's principle has 
proved valid under a very general metric (Kovner 
1990; Perlick 1990), we should be able to re-derive 
lens equation from our time delay result. Fermat's 
principle asserts that a light path is true if and 
only if its arrival time; or, what amounts to the 
same, its time delay is stationary with respect to 
variation of the turning point. In other words, 
d{/\t)/dgg = 0. With this sufficient and necessary 
condition for Fermat's principle, we can obtain the 
lens equation in the MOND regime. 

= - DLsy^l^. (60) 

V'^ is composed of three parts; and when g is very 
large 



VV'GJi = 4 



V^,^ = 27r 



Bo 

l-|-fc/4S7r 



1/2 



GNmg 
1 + ^nE/k 



(61) 

(62) 

(63) 



Obviously the addition of the last three equations 
equals A(f, given by Eq. (33), hence Eq. (60) is 
identical to the lens equations given by Eq. (34). 

5. Discussion 

Since the appearance of Bekenstein's relativistic 
MOND gravitation, the long-time incompleteness 
of Milgrom's modified Newtonian dynamics seems 
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to be filled up. Now we can investigate the GL 
phenomenon, such as the deflection angle or time 
delay, in the MOND paradigm, which could not 
be done before. 

In this paper,we have investigated the GL phe- 
nomenon under the approximations and presup- 
positions of 

(i) The GL lens is assumed to be static, spher- 
ically symmetric and following the thin lens for- 
malism, 

(m) The motions of light rays arc described in 
the framework of a Schwarzschild lens (i.e. a point 
mass model), 

{in) The revised physical metric is obtained by 
adding a positive scalar field into the potential of 
the standard Schwarzschild metric in symmetric 
coordinates. 

Under these presuppositions, we find that when 
9 is larger than 9^ (= tJOl] the MOND length 
scale), the reduced deflection angle a will ap- 
proach to a constant for a given mass. This spe- 
cial prediction is a feature of GL in TeVe S, which 
is similar to that obtained by Mortlock & Turner 
(2001) with an intuitional approach. This is no 
surprise because just like in GR, the deflection of 
photons is simply twice the deflection of a mas- 
sive particle with the speed of light in TeVeS, 
and thus Mortlock & Turner (2001) started from 
the correct premise. However, there is still some- 
thing that was unknown before the appearance of 
TeVe S. For a static spherically symmetric space- 
time, the case ^ 1, which yields a relation of 
Eq. (29), is valid only when |V$| <C (47r/fc)2oo. 
This condition is valid in the MOND regime when 
|V$| goes up to a couple of orders above ao, or 
equivalently go is around one order of magnitude 
below tq. (Bekenstein 2004) 

We should address that the criteria of mass dis- 
crepancy for GL effects in TeVe S consists with 
that of stellar dynamics. In other words, only 
when ^ kr^/A-K, the "missing mass" shall ap- 
pear. This corresponds to the demarcation of 
the high surface brightness (HSB) and low sur- 
face brightness LSB galaxies from the dynamical 
analysis (Sanders & McGaugh 2002). 

When we apply the deflection angle law to mag- 
nification, we find that in TeVe S the difference 
in the magnifications of the two images in the 
point mass model depends on the lens mass and 



source positions, and is always larger than one. 
This differs from traditional gravitational lensing, 
which says that the difference must always be 
one. Tens of thousands of the multiple images 
lensings found by the Sloan Digital Sky Survey 
(SDSS)(Stoughton et al. 2002) might be applied 
to check this prediction. For microlcnsing, light 
curves in TeVe S at deep MOND regime also differ 
from that in GR. To observe the discrepancy, the 
sources have to be located about Zs>'\- (Mortlock 
& Turner 2001). 

Concerning time delay, the result is even more 
exotic. Add an arbitrary positive scalar field 
into the primary Schwarzschild metric Bekenstein 
(2004) the resultant contributions of the scalar 
field will reduce rather than enhance the poten- 
tial time delay. Unfortunately, this phenomenon 
is unmeasurable in GL systems. What we can de- 
termine is only the time delay between two images 
produced by the same source. Therefore, the op- 
posite feature of the time delay in TeVe S, which is 
the same for all images, will be canceled out, and 
only those parts contributed from the defiection 
potential can be observed. 

Even though the opposite effect on the poten- 
tial time delay due to the scalar field can not be 
measured in GL system, the time delay between 
two images offers another constraint on the needed 
mass in GR and TeVe S, which usually differs from 
that given by the deflection angle. Actually the 
mass ratio obtained from these two approaches 
will be the same, otherwise the mass ratio given 
by time delay would always be smaller than that 
by deflection angle. In one word, the MOND and 
CDM paradigm are not mTitTially alternatives in a 
GL system, when we consider defiection angle and 
time delay at the same time. 

The first time delay for the gravitational lens 
Q0957-I-561 was measured in 1984 (Florentin- 
Nielsen 1984), since then more than 11 time de- 
lay lenses have been found, including 7 systems 
with a good quality of the astrometric data and 
2 systems with serious problems (Kochanek & 
Schechter 2003). It would be very interesting 
if those debatable systems can be explained in 
the framework of TeVe S. However, we emphasize 
that the measurement of time delay is controver- 
sial. Actually it took almost 20 years of debate 
between a short delay and a long delay on the 
first observed time delay source (Fetters, Levine 
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& Wambsganss 2001). 

It also worth noting that our conclusion is based 
on the Bekenstein's approach to investigate a grav- 
itational lens. He assumes a priori that the po- 
tential dominating in the massive particle dynam- 
ics agrees with the potential influences on the GL 
system. In other words, the temporal and spa- 
tial components of the metric respectively have a 
form -(1 + 2<I>) and (1 - 2<1?), whence gugu ~ 1- 
However, there is still a chance that these two po- 
tentials are not the same and even that gugu 7^ 1 
under the framework of TeVe S. If so, the conclu- 
sion of this paper may be different (Edery 1999; 
Bckenstein et al. 2000). 

Since we can observe the influence of gravity 
on gravitational lens phenomena at distances fur- 
ther away than we can for the dynamics of massive 
particles, GL systems offer a good method to dis- 
tinguish between GR and TeVeS. Although we 
need more than the point mass model to fit astro- 
metric data, this simplified model has in principal 
told us the exotic predictions of the deflection an- 
gle and time delay. It would be very interesting if 
we could develop a more realistic model to fit the 
known data in GL systems. 
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